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The study and characterization of the diversity of spatiotemporal patterns generated when a 
rectangular layer of fluid is locally heated beneath its free surface is presented. We focus on the 
instability of a stationary cellular pattern of wave number ks which undergoes a globally subcritical 
transition to traveling waves by parity-breaking symmetry. The experimental results show how 
the emerging traveling mode {2/^ks) switches on a resonant triad (A;^, /cs/2, 2ks/3) within the 
cellular pattern yielding a "mixed" pattern. The nature of this transition is described quantitatively 
in terms of the evolution of the fundamental modes by complex demodulation techniques. The 
Benard-Marangoni convection accounts for the different dynamics depending on the depth of the 
fluid layer and on the vertical temperature difference. The existence of a hysteresis cycle has been 
evaluated quantitatively. When the bifurcation to traveling waves is measured in the vicinity of 
the codimension-2 bifurcation point, we measure a decrease of the subcritical interval in which 
the traveling mode becomes unstable. From the traveling wave state the system undergoes a new 
global secondary bifurcation to an alternating pattern which doubles the wavelength {kg/ 2) of the 
primary cellular pattern, this result compares well with theoretical predictions [P. Coullet and G. 
loss, Phys. Rev. Lett. 64, 866 (1990)]. In this cascade of bifurcations towards a defect dynamics, 
bistability due to the subcritical behavior of our system is the reason for the coexistence of two 
different modulated patterns connected by a front. These fronts are stationary for a finite interval 
of the control parameters. 

PACS numbers: 47.55.pb, 45.70.Qj, 47.20.Ky, 47.20.Lz 



I. INTRODUCTION 

The experimental and theoretical study of the dynam- 
ics of ID cellular patterns is still a subject of great in- 
terest, specially when this research is devoted to spa- 
tially extended systems. The simplicity of such sys- 
tems, combined with the theoretical approach by nonlin- 
ear dynamics and critical phenomena, allows us to dis- 
cover universal mechanisms at the characteristic length 
scales without being necessary to cope with the com- 
plex microscopic problem. These ID dynamics for spa- 
tially extended dissipative systems have been reported 
for a large variety of experiments where the bifurcation 
scenario from a homogeneous pattern can be a primary 
bifurcation either to a stationary cellular pattern or to 
a traveling wave pattern. As the control parameter is 
increased the system may go through a secondary bifur- 
cation. Some of these experiments are: the printer's in- 
stability [1, 2j, circular liquid column arrays O [H |5], 
directional solidification [6 , directional viscous finger- 
ing 0, the Taylor-Dean system [8], electroconvection in 
liquid crystals with rectangular cell [9l [10] , lateral heat- 
ing in a rectangular cell (TTJ [121 [131 [TH [15] , lateral heat- 
ing in an annular cell [16], ID heating in a rectangular 
cell [ni[Tl|ll|20l[2ll|2a[2a[24], numerical simulations 
with the stabilized Kuramoto-Sivashinsky equation [25] . 

The aim of this paper is the study of a dissipative ex- 
tended system, a buoyant-thermo-capillary driven flow, 
which undergoes a cascade of secondary bifurcations 
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from a stationary cellular state towards a defect dynam- 
ics. Certain symmetries are broken in each bifurcation. 
Specifically, we study experimentally the mode competi- 
tion at the threshold of a secondary instability to travel- 
ing waves. This convective approach allows us to carry 
out experimental measurements where the relevant mag- 
nitudes can be easily measured, with the cost of relatively 
large time-scales required to achieve asymptotic states. 

Our system consists of a thin rectangular fluid layer 
of silicone oil opened to the atmosphere with an inho- 
mogeneous heating at the center and along a line placed 
at the bottom. This experimental set-up is intended to 
approach an "ideal" quasi- ID thermoconvective system. 
The fluid dynamics is in the context of the Oberbeck- 
Boussinesq approximation for the applied range of tem- 
peratures. The physically relevant control parameters are 
the vertical temperature difference AT^ and the depth of 
the fluid layer d. 

Without threshold, as soon as AT^ ^ 0, a homoge- 
neous convective pattern (PC) appears: the fluid over 
the heater line (HL) raises by the effect of the negative 
vertical temperature gradient and it returns down near 
the cooling walls. Therefore, a circulating flow at both 
sides of HL develops the structure of two counterotating 
rolls whose axis are parallel to the HL. As we increase 
the control parameter AT^ from PC, this pattern will 
invariably remain for the subsequent instabilities. For 
asymptotic time responses we can distinguish three fun- 
damental patterns: a stationary cellular pattern (ST) 
with wave number /c^, a traveling pattern (TW) consist- 
ing of one propagative wave with wavenumber 2ks/?> and 
an alternating pattern (ALT) consisting of two traveling 
counterpropagative waves with wavenumbers kg/ 2 plus 
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the stationary one (a resonant triad). Similar patterns 
have been analyzed in a previous experimental setup [17^, 
but with a less developed system. In this work, accuracy 
has been improved by means of new optical and data 
acquisition systems. 

In 1990, Coullet and loss [2F (CI) classified the global 
character of secondary bifurcations in ID systems using 
a Bloch-Floquet analysis; and predicted ten generic sec- 
ondary instabilities. In this theory, there is a coupling 
between the phase of the basic pattern and the amplitude 
of the unstable modes associated to the secondary bifur- 
cation. On the other hand, the phase-amplitude coupled 
equations obtained with less restrictive symmetry argu- 
ments by L. Gil [27l[28] are suitable for experiments with 
analogous dynamics, but six real and fifteen complex co- 
efficients are necessary to be determined experimentally. 

We report the instability of a "mixed" pattern due to 
the close position of the new emerging mode of the TW. 
This mixed pattern consists of localized ALT domains 
with ffuctuating boundaries over ST. These results shed 
new light on the role of two close competing oscillatory 
modes in the onset of a subcritical transition to TW. 
These are the emerging traveling mode and the nearest 
oscillatory one with wavenumbers 2ks/?> and kg/ 2^ re- 
spectively. This relationship dos not fit into the classifi- 
cation made by CI, and this fact could be the consequence 
of (z) subcriticality or (zz) the dynamical implication of 
the other two spatial directions. 

In the bistable regimes for two subsequent subcriti- 
cal bifurcations to TW and ALT, we report existence 
of normal fronts [29 that connect two stable states of 
the oscillatory type. These fronts are stationary and ac- 
cordingly to the results reported here, they are pinned 
to both modulated patterns for a finite control parame- 
ters range. A similar front dynamics linked to subcriti- 
cality has already been modelized for binary fiuids [30 , 
however in this case pinning is produced between a ho- 
mogeneus and a traveling wave pattern. In transients, 
when we place our system below and close to the tran- 
sition point, we report the study of the so called Fisher- 
Kolmogorov-Petrovsky-Piskunov (FKPP) fronts [31 that 
connect stable states with unstable states. These fronts 
provide interesting features of the dynamics such as the 
convective vs. absolute character of an instability. 

From the theoretical point of view subcritical transi- 
tions are predicted in superconductors and in liquid crys- 
tals [32]. Despite subcriticality has been widely studied 
experimentally in other systems, some examples (not pre- 
tending to be an exhaustive list) are: pearling of tubular 
membranes [33], the Freedericksz transition [34] and the 
transition to traveling waves in binary fiuids [35], here 
we report quantitative measurements for a new subcrit- 
ical transition in a quasi- ID convective system from a 
cellular pattern to traveling waves and further beyond. 



II. EXPERIMENTAL SETUP 

The fiuid layer of depth d is placed in a narrow con- 
vective channel x Ly {L^ = 470 mm, Ly = 60 mm). 
As shown in Fig. [l^ this layer lies over a fiat surface, 
a mirror of thickness 3 mm. The bottom of the mir- 
ror is in contact with a heater rail of 1 mm (thickness) 
which provides the heating line in the x direction, HL. 
The heater rail belongs to an aluminum block with an 
inner closed water circuit thermoregulated by a heater 
bath at T^. The temperature measured at HL on the 
mirror surface differs from the one given by the heater 
bath probe in -3.00 °C for d = 4 mm and Th = 48 °C. 
We limit the heat transfer in the y direction by isolat- 
ing laterally the heater rail with Plexiglas. Therefore in 
the transversal direction of the cell, one may approxi- 
mately achieve a Gaussian temperature profile smoothed 
nearby the central heater because of the thermal diffusiv- 
ity of the mirror. Regarding other experimental setups 
like the resistive wire experiments (20 jam - 50 /im di- 
ameter) [21] [22] [23^ where there is a coupling between 
the circulation ffow and the thermal stability of the wire, 
in our system because of the great thermal inertia cou- 
plings between the convection in the bulk and the heater 
are not expected. The lateral walls are two aluminum 
blocks (coolers) whose temperatures are kept constant at 
Tc = 20 ± 0.1 °C by means of a secondary water cir- 
culation. This temperature is given by the cooler bath 
probe and it is the same as the one measured at the top 
of the aluminum blocks. The shorter boundary walls at 
the opposite extremes on x are made of Plexiglas. All 
these experimental parts are contained in a Delrin block. 
The effective dimensions of the cell are 450 mm x 60 mm 
X 15 mm {Lx x Ly x Lz)^ although the optical setup al- 
lows us to visualize only a centered area from an upper 
surface of 310 mm x 60 mm. 

The geometric aspect ratios for d = 4 mm are: F^^ = 
L^/d = 112.5, F^ = Ly/d = 15. The ratio T^/Ty = 7.5 
allows us to classify the system as weakly confined in x. 

The fiuid used is a silicone oil (brand name: Dow Corn- 
ing 200) with a nominal viscosity of 5 cSt (see properties 
on Table [ij). This fiuid is transparent to visible light and 
therefore it is suitable for the shadowgraph technique. It 
has a low volatility, the depth of the fiuid layer diminishes 
with a rate of 0.01 mm/24 hours over an area of 270 cm^ 
at 20 °C, so we can work (in the range of temperatures 
applied) with the layer opened to the atmosphere. Un- 
der this condition, we assume that the depth, measured 
with a precision of 0.01 mm with a micrometric screw, 
remains constant along each measurement. The temper- 
ature outside the cell is controlled at = 20 ± 0.1 °C 
for a sequence of measurements. 

From the corresponding physical properties (Table l| 
there are two significant adimensional numbers: the 
Prandtl number Pr = v / n ^ 75 (= compares viscous 
diffusivity vs. thermal diffusivity) and the dynamic Bond 
number Bojj = R/M ^ 1.05 {d = 3 mm) (= compares 
thermogravitatory effects vs. tliermocapilarity effects). 
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FIG. 1: Sketch of the experimental setup, (a) Rectangular 
cell: cross section (up) and top view (bottom) of the rectan- 
gular vessel. The scaling in both views is different; (b) Optical 
setup. 



TABLE I: Physical properties of the 5 cSt silicone oil. 



Surface tension (25° C) 
Thermal conductivity (50° C) 
Thermal diffusivity (A/pCp) 
Kinematic viscosity (25° C) 
Density (25° C) 
Refractive index (25° C) 
Linear expansion coefficient 
Surface tension/temperature 



a ■■ 

A: 



19.7 mN m"^ 
0.117 W m-^ K-^ 
K = 6.68 10"^ m^ s"^ 

= 5 cSt (5 10"^ m^ s" 
p 913 kg m"^ 
n = 1.3960 
a = 0.00105 K-^ 
1^ = -8 10-^ N m-^ K- 



As well ai d = 3 mm and ATy = 30 K, we can deter- 
mine two characteristic time scales: the viscosity time 
scale Tjy = 2 s and the thermal diffusivity time scale = 
135 s, thus the dynamics is driven by the temperature 
field. 

The shadowgraphic flow-visualization system is 
sketched in Fig. In order to obtain quantitative 

results for a sequence of measurements, the light 
intensity from an incandescent light bulb (white and 
incoherent light) is controlled with a power supply. An 
incoming parallel light is sent through the convection 
pattern. The light beam crossing the fluid layer suffers 
different deflections depending on the refractive index 



gradients (i.e. temperature gradients). The analysis of 
the spatiotemporal periodicity of the modulation of the 
light beam, that has crossed the fluid layer, provides 
information about the temperature field of a certain 
pattern. In a similar experimental cell the deflection 
induced by traveling waves was below 1 /im for a depth 
of 3 mm [17 , therefore the information analyzed comes 
from the convective flow in the bulk. Once it is reflected 
back at the mirrored bottom the output beam is pro- 
jected into a screen. The shadowgraph method allows us 
to observe the thermal gradients in the bulk of the fluid 
layer as a modulation of the light intensity on the screen. 
The screen image is recorded with a CCD video camera 
with a resolution of 570 px x 485 px (px = pixels), 
keeping constant the focal length, the aperture and 
the gain. The shadowgraphy image is an instantaneous 
image of the screen (Fig. [2|. There, we can observe a 
central dark line (HL) and a bright modulation in both 
sides of HL by symmetry due to the effect of thermal 
lenses. Periodically, a line over the bright modulation 
next to HL is recorded to obtain a spatiotemporal dia- 
gram 6'(x,t), and another line placed perpendicular to 
HL is also recorded to obtain a spatiotemporal diagram 
S{y^t). With the developed image analysis software, the 
signal S{y,t) is averaged over a temporal sequence to 
determine a medium contrasted position in y to place 
the acquisition line that records the spatiotemporal 
diagram S{x^t) (with an acquisition rate of 1 s~^). The 
signal S{x,t) is processed via a bidimensional Fourier 
transform and complex demodulation techniques in 
order to determine the amplitudes Aj , the wave numbers 
kj and the frequencies ujj of each fundamental Fourier 
mode M(kj^uOj) [36]. S{x,t) can also be filtered by 
discriminating the normalized amplitude below a critical 
value /i. 



The control parameters are the depth of fluid layer 
d and the vertical temperature difference, ATy = — 
Ta. Under these conditions, we have obtained the smaller 
threshold for the primary bifurcation at d = 4.5 mm. To 
build the stability diagram we have explored depths d = 
2.5, 3, 3.5, 4, 4.5 and 5 mm. To locate the different 
asymptotic regimes we have worked with step |AT^| = 
2 K, except at d = 2.5 with a bigger step {\ATh\ = 5 
K) and some isolated measurements at d = 3.5 looking 
for TW. Meanwhile boundaries have been defined with a 
maximum step of \ATh\ = 1 K. To determine the nature 
of the bifurcations at d = 4 and 4.5 mm, we increase and 
decrease ATy by steps \ATh\ = 1, 0.5 and 0.3 K, allowing 
the system to achieve the asymptotic state for 1 hour in 
the ascending sequences, and 3 hours in the descending 
sequences nearby the threshold. Experimentally, we have 
determined that the relaxing time for a step of AT^ = 1 
K is approximately Ih and the relaxing time for a step 
of ATh = -1 K is approximately 2 hours and 30 min. 
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FIG. 2: Shadowgraphy images for three different regimes, 
below the dynamics of hotspots over HL is sketched (a zoom 
view of the selected area above): (a) stationary state ST for 
d = 3 mm; (b) traveling waves TW for d = 4 mm; (c) al- 
ternating waves ALT for d = 4 mm. The arrows between 
hotspots indicate the direction of their movement. 



III. RESULTS 
A. Dynamics of the hotspots 

The effect of thermal lenses produces on the screen a 
convergence of the emerging rays corresponding to the 
bright profiles (see Fig.[2|, as well as a divergence for the 
ascending convective flow over HL which corresponds to 
a distribution of aligned dark spots or hotspots. In Fig. |2] 
we show the shadowgraphy images belonging to the pre- 
viously introduced regimes: ST (stationary), TW (trav- 
eling wave) and ALT (alternating). Below the shadowg- 
raphy images we have sketched how the dynamics of the 
hotspots is observed over HL and the space symmetries 
of the contours of maximum brightness. 

If we take the wavelength of the basic cellular pattern 
ST, Ag, as the space occupied by a hotspot, we have 
got roughly 75 oscillators which are differently coupled, 
regarding the type of pattern, in a dissipative system. 
At the same time, a driving force is applied into them 
when we increase ATy. The ratio between wavelengths 
of different patterns is conserved for different depths (i.e. 
see Figjsk for d = 4 mm). The frequencies are shown 
in Fig. |8p, and the phase velocities can be obtained as 

= Lo/k. 

For the secondary instabilities that are going to be 
quantified in the following sections we show in Fig. [3] the 
characteristic spatiotemporal diagrams 5'(x, t) and their 
respective bidimensional Fourier spectra {k^co). The fun- 
damental modes for each pattern have been pointed out 
as Ms (the stationary mode) and My± (the traveling 
modes to the right (+) and left (— )). The dynamics 
of the hotspots allow us to describe the different regimes: 



(i) The continuous spatial translation symmetry of the 
pattern PC is broken into a steady cellular pattern 
ST (Fig. [3^). In this pattern the hotspots keep 
their positions fixed at a distance Ag ~ 2d (i.e. at 
d = 4mm see Fig. |8^) of their first neighbors, so 
we have a discrete translation symmetry on space 
{x) and a continuous one on time (Fig. [2|i). The 
corresponding fundamental mode is Ms{ks^O). 

(ii) The transition from ST to TW breaks the parity 
symmetry {x —x) and the temporal continuous 
translation symmetry. In the TW regime (Fig.jsjD), 
the new discrete translation symmetry on space 
is given by a wavelength Xtw ^ 3/2 A5 = A35/2 
(i.e. at d = 4mm see Fig. [8^). This result differs 
from the doubling period obtained in a previous 
work [17] because the precision for both, the experi- 
mental techniques and the analytic data processing, 
has been remarkably improved. The hotspots de- 
cide to travel at the same time towards a privileged 
direction selected by the system at a frequency 
ootw^ which varies linearly with ATy (Fig. [s})). 
The mean period is approximately 25 s. The ve- 
locity at which hotspots travel is the phase veloc- 
ity V(/y = A35/2 ^TW^/27r, because we have measured 
that the linear group velocity verifies Vg = duj/dk ^ 
0, as it will be shown further on. So the phase 
velocity increases linearly with an increasing con- 
trol parameter ATy as we go beyond the thresh- 
old of the secondary bifurcation to TW. The corre- 
sponding fundamental mode is (either right or left): 

My±{k2s/3,^^TW)' 

(iii) In the ALT pattern (Fig. ^) the hotspots are oscil- 
lating in counterphase with a new discrete spatial 
translation symmetry given by Xalt ^ 2As = A2S 
(i.e. ai d = 4mm see Fig. Isk) in agreement with 
previous experimental work"lT]. So hotspots be- 
have as oscillators that by pairs interchange their 
positions to recover their initial positions after a 
cycle of frequency coalt- This frequency is almost 
constant as we increase ATy (Fig. jsj)). The mean 
period is about 23 s. So, for an ascending sequence 
the phase velocity keeps constant as we go beyond 
the threshold of the secondary bifurcation to ALT. 
This pattern is strongly nonlinear as the result of 
the competition between the fundamental modes of 
a resonant triad: Ms{ks , 0) and My± {ks/2 , ^^alt)- 
The stationary mode Ms{ks^ 0) characteristic of the 
cellular pattern ST has been a sort of ^'restored\ 

In this array of coupled oscillators there may exist 
localized discontinuities that correspond to boundaries. 
These boundaries or fronts break the dynamics into two 
different groups of oscillators, each one of them with a 
great spatiotemporal coherence. On the spatiotempo- 
ral diagram S{x^t) these regions are called localized do- 
mains. Such is the case of the existence of one or various 
localized domains in the ALT pattern coexisting with the 
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FIG. 3: On the left the most representative spatiotemporal 
diagrams S{x,t), and on the right the contour plot of the cor- 
responding Fourier spectrum {k,ou) (the fundamental modes 
Ms (stationary) and Mv± (traveling) have been pointed out, 
the central spot corresponds to (/c = 0, a; = 0)): (a) ST (d = 
3 mm, ATy = 16.7 K) ; (b) TW (d = 4 mm, ATy = 18.6 
K); (c) ALT (d = 4 mm, ATy = 25 K); (d) mixed pattern 
ST/ ALT (d = 3 mm, ATy = 30 K), the boundary of the lo- 
calized pattern ALT has been highlighted (filtering with /x = 
0.37). 

ST pattern as a "mixed" state (Fig. [sji). This domains 
usually collapse, they have fluctuating boundaries with 
variable spatiotemporal areas. This pattern from now on 
is designated as ST/ ALT. 

B. Primary and secondary instabilities 

A smooth increase of the control parameter ATy can 
lead from the homogeneous pattern PC to the onset of 
two different kinds of primary bifurcations depending on 
the depth of the layer d (see the stability diagram in 
Fig:[4|: 

(i) For d < 3.5 mm, as soon as ATy ^ the ST pattern 



becomes unstable. For this pattern we find a higher 
threshold as we increase the depth d from 3.5 mm 
to 4.5 mm. 

(ii) For (i > 4.5 mm, the TW pattern becomes unsta- 
ble from the primary roll PC. This bifurcation is 
thought to be supercritical concerning the experi- 
mental fact that the TW pattern appears simulta- 
neously along the cell from an homogeneous pattern 
at d = 5 mm. 

For the explored depth range d = 2.5-5 mm, there are 
the following pathways to secondary bifurcations (see the 
stability diagram in Fig:[4|: 

(i) Thick layers: For 4.5 < <i < 5 mm and further 
from the threshold of the primary bifurcation to 
TW, the system undergoes a new bifurcation to 
the ALT pattern. The relative amplitudes verify 
\Ay_\ ^ \Ay^\ and ^ |A,±|/2. 

(ii) Intermediate layers: For 4 < d < 4.5 mm, 

the system becomes unstable from ST to TW via 
a bistable regime. This bistable regime consists 
of a mixed pattern, ST/ ALT, coexisting with the 
new TW pattern. The front connecting these two 
regimes (ST/ ALT and TW) is stationary for the 
front velocity is Vp ^ 0.01 mm/s. This front ve- 
locity is determined by measuring the slope of the 
front defined by the new mode, which becomes un- 
stable, at the spatiotemporal diagram. In order to 
study this bistability we focus on this second insta- 
bility to TW for two intermediate depths: 4 and 
4.5 mm. Beyond the TW pattern we report a new 
secondary bifurcation towards the ALT pattern at 
d = 4 mm. 

(iii) Thin layers: For d < 3 mm, from the basic flow 
pattern ST, a cascade of secondary bifurcations 
takes places towards a pattern of localized domains 
over the basic pattern ST. These results will be 
reported elsewhere. The first instability from the 
basic pattern ST is the mixed ST/ ALT pattern 
(Fig. [sji). The quantitative analysis shows that 
\Ay±\ keeps constant inside these domains along a 
sequence for increasing steps. 

For thinner layers {d < 2.5 mm) and high temperatures 
{ATy > 55 K) the heating looses its linear localization 
and hence it turns into an homogeneous 2D heating. Un- 
der these conditions, the system gives birth to an hexag- 
onal pattern arranged symmetrically with respect to HL, 
similar to the classical Benard-Marangoni convection. 

Beyond the threshold of the last secondary bifurcation, 
for each d and until the maximum value of ATy applied, 
the dynamics of the system still keeps the characteristic 
patterns, but with the presence of defects of the type 
of dislocations with topological charge +1. For transient 
regimes we have observed the existence of pulses, sources, 
sinks and localized drifting domains (see Fig.[5|. Also, in 
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transients we can find counterphase oscillating patterns 
(optical modes). For steady patterns when there is a 
domain boundary (or front) connecting the patterns TW 
and ST/ ALT, these oscillating patterns are present as 
well. As it is shown later on, oscillatory patterns are 
similar to the ALT pattern on account of the fundamental 
modes but with different amplitude rates. 

From now on, our aim is to study the subcritical nature 
of the system in the secondary bifurcation towards TW 
(with a backward study research looking for an hysteresis 
cycle) and towards ALT. 



C. Secondary bifurcation to traveling waves and 
further 



For intermediate layers, when the system approaches 
the threshold from the cellular pattern ST to the sec- 
ondary bifurcation to TW we define the reduced control 
parameter as e = ATy/AT^^ — 1. AT^^ is defined as 
the critical parameter at which the system becomes un- 
stable towards the new state, TW, for an ascending se- 
quence of measurements with the minimum step \ATh\. 
The amplitude of the corresponding unstable mode \Ay-\ 
grows abruptly at this point meanwhile the amplitudes 
I As I and decay to zero (i.e. see Fig. The 

subcritical parameter for the hysteresis cycle is given 
by Ec = AT^JAT^^ - 1 < 0, where AT^^ is the criti- 
cal parameter (corresponding to a descending sequence) 
at which the system returns to the original state, ST. 
Thus, the range Sc < s < is the subcritical interval. 
At d = 4 mm, for an ascending sequence with step 0.5 
K from the amplitudes diagram (Fig. [7^) we have mea- 
sured AT^^ = 17.7 K. To verify the existence of hysteresis 
we have worked with ascending and descending steps of 
0.3 K in the subcritical interval and we have obtained 
AT^^ = 16.4 K and Sc = -0.07. 

Close to the codimension-2 point (at d = 4.5 mm) de- 
scribing a cycle with step 0.3 K and with the same pro- 
cessing protocol we have determined Sc = -0.02. 

In a correlative ascending sequence for the secondary 
bifurcation to TW (Fig.[6|: 

- For Ec < e < the mixed pattern ST/ ALT is 
present. In the localized domains of the ALT pat- 
tern the competition between the traveling modes 

My^{ks/2,usT/ALT) and My_{ks/2, -cjst/alt)^ 
with similar amplitudes, produce a resonant 
nonlinear interaction with the stationary mode, 
M,(/c„0). 

- For e = the system bifurcates towards TW at the 
same time that coexists with the original ST/ ALT 
pattern, this is the bistable regime. In this regime 
a discontinuity of the amplitude of the traveling 
mode M^_(/c2s/3, —^tw) has been produced in the 
corresponding domain, the TW pattern. The front 
velocity is approximately zero, Vp ~ 0.01 mm/s. 



- For £ > the whole system has bifurcated to TW 
without the existence of boundaries, neither sources 
nor sinks, so M^_(/c2s/3, —^tw) is a global unsta- 
ble mode. 

For this secondary bifurcation to TW if we artificially 
excite (perturbing locally the free surface) the basic state 
ST just below -Sc at (d = 4 mm, ATy 15.5 K), we 
observe the birth of a TW domain (see Fig. [9^). This 
domain vanishes after a while because of the attenuation 
of the amplitude of the unstable traveling mode (Fig. [oJd) 
with a front velocity of Vp ~ 0.17 mm/s. This is what 
is expected for a negative growthrate when we are below 
and close to the threshold. For these control parameter 
values the TW pattern is said to be convectively unsta- 
ble with a negative group velocity: \vg\ ~ 0.12 mm/s, 
and then, the phase velocity verifies ~ -3.5 Vg. If we 
perturb the system above and close to the threshold, the 
group velocity is null. 

Further from threshold at d = 4 mm the basic state 
of TW bifurcates to the new state, ALT, via a dynam- 
ics of stationary fronts at the asymptotic states (see the 
spatiotemporal and filtered diagrams in Fig. 10i,b). For 
these fronts we have measured a front velocity of Vp ~ 
0.01 mm/s. We define a new threshold for an ascending 
sequence (with step IK) as e"^ = ATy/ ATy c — 1 = 
where at ATyc the ALT pattern appears for the first 
time {ATy = 21.5 K). From this point on and as far 
as we increase AT^, the front invades the TW pattern 
until at e'^ = 0.29 the whole system dynamics has tran- 
sitioned to the ALT state (Fig. 10 3). In this transition 



we have measured a finite jump for the amplitudes of 
the traveling modes at the threshold: \Ay±\{ALT) ^ 
3\Ay-\{TW) while the amplitude of the stationary mode 
fulfils: \As\{ALT) ^\Ay±\/2{ALT). 

For higher |AT(;|, the turbulent ALT pattern has a dy- 
namics of defects (dislocations) (Fig. [lT^,b). Buoyancy 
plumes can be observed in the shadowgraphy image, they 
are detached along the perpendicular axis to HL. On the 
phase gradient in Fig. pTj : we can observe how these de- 
fects are spontaneously generated, adjusting locally the 
wavenumber of the pattern. 

Because of the subcritical behavior of the system in the 
secondary bifurcation to TW, nonlinear interactions be- 
tween the original (ST/ ALT) and the new (TW) patterns 
in the bistable regime (Fig. 12 1) are expected. Indeed 
the boundary is likely to show this nonlinear behavior 
by resonant interactions between the competing modes. 
In consequence, we focus not on the evolution of each of 
the coexisting domains, but on the overlapping between 
the original mode Mg (representing the mixed ST/ ALT 
pattern) and the new mode My_ (representing the new 
pattern TW) at the front. For this analysis the signal 
S{x^t) has been filtered with /i = 0.37. The average 
width of the front is defined as a subcritical length, Lsub^ 
where the mode Mg is interacting nonlinearly with the 
others {My±). 

To explain the resonance of the amplitude \As\ as 
a nonlinear coupling between the triad (M5, My- and 
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d(mm) 



FIG. 4: Stability diagram (d, ATy) built from the spatiotemporal diagrams S{x,t). PC and ST are stationary patterns, TW 
and ALT are oscillatory patterns and ST/ ALT is the mixed pattern. ST/ ALT + TW defines the subcritical region where a 
stationary front connects ST/ALT and TW. Double arrows correspond to the stepped sequences of measurements reported. 




FIG. 5: Spatiotemporal diagram of TW = 4 mm, ATy = 
20 K) with localized drifting domains in the ALT pattern in 
a transient (for an upward jump with step ATh = 4 K) and 
the presence of a source of TW in the left hand side. 



M^+) along Lsub^ we study the profile of the modulus of 
the amplitudes of the fundamental modes in the direction 
perpendicular to the front, so we can check the amplitude 
gain of \As \ at t = 131 s (see Fig, and the minimum 
value of 1 at t = 489 s (see Fig. ^) . The mode 
Ms is coupled with the propagative modes in the bound- 
ary where the unstable mode My- penetrates beyond the 
stationary front, and therefore in this region the modulus 
1 will be greater. Filtering with ja = 0.37 and for t = 
131 s, in the region Lgy^j^^ we obtain the following results 
|As|/|A^_| ^ 2 and |A^_|/|A^+| ^ 6. In the bound- 
ary the kinematic resonant conditions [37] for the wave 
numbers: 0.7774 mm~^ ^ 0.4438 mm~^ +0.3135 mm~^; 



and for the frequencies: 0.0030 s-^ ^ -0.2102 s"^ + 
0.2210 s-i are fulfilled. 

For each one of these data the corresponding error from 
the analytic processing is negligible considering the ex- 
perimental error due to temperature inhomogeneities. 

The corresponding oscillatory pattern that shapes the 
boundary contour is very similar to the previously men- 
tioned for the optical modes. This pattern differs from 
the ALT uniquely on account of the different ratios be- 
tween amplitudes. 

We can interpret the growth of the modulus \As\ 
(Fig. as a positive energy transference from the 

propagative mode My- (TW) to the stationary one Ms 
(ST/ ALT). At t = 489 s the amplitude of Ms has been 
attenuated regarding the previous result at t = 131 s, 
meanwhile for the propagative modes the ratio of ampli- 
tudes accomplishes: |74^_|/|A^+| ^ 1 (see Fig, p^). 

Concerning the correlation length for the subcritical in- 
stability to TW, we may associate the subcritical length 
Lsub with an attenuation length of the new pattern (TW) 
inside the original one (ST/ ALT). Nevertheless, the evo- 
lution of the averaged front width, as a correlation length 
parameter, has not revealed any evidence at the onset of 
this transition. Fluctuations of the amplitude \As\ are 
caused by the nonlinear resonance at the front, so this 
fact has to be considered to implement the most suitable 
strategy for the data analysis in bistable regimes. 

IV. DISCUSSION AND CONCLUSIONS 
A. An extended ID convective system 

The two counterotating rolls (PC) in this experiment 
remain for any AT^, so we might compare the basic cel- 
lular pattern ST (which is a transverse cellular pattern 
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FIG. 6: Spatiotemporal diagrams corresponding to the secondary bifurcation to TW at d = 4 mm (for an ascending sequence): 
(a) the mixed ST/ ALT pattern belongs to the subcritical interval, the boundary of the ALT domain has been highlighted 
(filtering with /x = 0.5); (b) the bistable pattern: ST/ALT and TW at the threshold, these patterns are connected by a 
stationary front; (c) the TW pattern above the threshold. 




FIG. 7: (a) Bifurcation diagram at d = 4 mm for ascending steps of 0.5 K. The continuous line is a guide to the eye. The 
discontinuous arrow show the upward jump of the amplitude mode Mv- at ATv = 17.7 K (s = 0) for the secondary bifurcation 
to TW. Far from this threshold at ATy ~ 19.2 K the amplitude of the mode My- decays in favor of the counterpropagative 
mode Mv-\- representing the advance of the front which splits the patterns TW and ALT. (b) Sketch of the mechanism of 
activation of resonant triad in the ALT domains (mixed ST/ ALT pattern) by the emerging TW modes, ujsr is the frequency 
at the subcritical interval {—Sc < £ < 0). 



regarding the wave vector of the primary roll) with other 
systems with similar behavior such as the lateral heat- 
ing experiment in a rectangular cell with Pr = 10 [11 . 
In this experiment, the fluid dynamics for deep layers 
corresponds to stationary rolls with a wavevector per- 
pendicular to the temperature gradient, this pattern be- 
longs to the type of stationary rolls predicted by Smith & 
Davis [38]. These stationary rolls become unstable over 
the underlying "unique" roll which invades the whole cell. 
For larger horizontal temperature gradients the system 
becomes unstable via an oscillatory pattern of optical 
modes. This pattern has a strong affinity with the reso- 
nant triad of the ALT regime where the counterpropating 
modes My± have restored the stationary mode Mg. 
For each pattern, the dynamics of the central 31 cm of 



the cell is equivalent to the dynamics of an array of 56 os- 
cillators (from a total of 75) interacting with their neigh- 
bors. We are studying the collective behavior of a set of 
macroscopic particles (hotspots), a ID array of nonlinear 
coupled oscillators that become unstable. Experimen- 
tally we are performing "coarse-graining" when we asso- 
ciate to the hotspots on HL the microscopic variables of 
the convective system like the thermal fluctuations. Be- 
sides, as the physical aspect ratios for the cellular pattern 
at d = 3 and 4 mm are respectively Tx = L^/Xg = 75 
and 59, we may consider our system to be extended. 

According to the experimental results the system may 
suffer two different bifurcations to TW: subcritical and 
supercritical, depending on the basic state and thus, on d. 
In consequence, from the homogeneous state PC the sys- 
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FIG. 8: At d = 4 mm evolution vs. ATv (with step ATh = 0.5 K) of: (a) the dimensionless wavenumbers (horizontal continuous 
lines correspond to the averaged wavenumbers for the same patterns at 3 mm and 4 mm) and (b) the frequencies {tu = 3.2 s) 
for each involved fundamental mode: Ms (stationary mode), Mv± (right-left traveling modes). Vertical dashed lines separate 
different patterns. 
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FIG. 9: (a) Spatiotemporal diagram with a localized mechanical perturbation close to the threshold to TW at (d = 4 mm, 
ATv = 15.5 K). The boundary of the transient TW domain (filtered with /i = 0.37) has been sketched; (b) Attenuation of the 
module of the traveling mode Mv- inside the transient domain. 




FIG. 10: Spatiotemporal diagrams corresponding to a secondary bifurcation to ALT above the threshold: (a) There is a front 
connecting the TW and ALT patterns at e'^ = 0.11; (b) Contrasted profile of the front (The left traveling mode has been filtered 
with /i = 0.37) at e"^ = 0.15; (c) the ALT pattern has invaded the whole cell far from the threshold at e'^ = 0.38. 
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FIG. 11: ALT pattern with the presence of defects at d = 5 
mm and ATy = 62 K: (a) The shadowgraphy image on the 
screen; (b) the spatiotemporal diagram; (c) phase gradient 
of the stationary mode Ms, the discontinuities correspond to 
dislocations. 



tern undergoes an oscillatory supercritical Hopf bifurca- 
tion that can be theoretically described by the cubic com- 
plex Ginzburg-Landau equation. Experimentally, this 
kind of supercritical bifurcation has been found in many 
other extended systems |lTl [24] and references therein. 
Meanwhile the system is behaving subcritically for d < 
4.5 mm. 

From these results we may infer that as the dynamic 
Bond number is Bod oc cP, when Bod decreases (as we 
move to the left in the stability diagram of Fig. |4| a sub- 
critical nature of the secondary bifurcations arises. So in 
this experiment subcriticality may be linked to the fact 
that thermocapillary effects are becoming more impor- 
tant the thinner the layer is. Subcriticality via a Hopf 
bifurcation to TW has been reported experimentally for 
example in binary mixture of fluids [39]. 

Our results for the secondary transition to TW show 
that Xtw/^s ~ 3/2. Under the light of this result we 
leave an open question to derive a suitable model as the 
one developed by L.Gil ^3 ^ system. This model, 
for a Floquet exponent q = 1/2, differs from CI because 
the constraint over equal amplitudes of the right-left en- 
velope functions (propagating modes) is relaxed. Hence 
this model can describe mixed patterns and localized 
coherent structures from a secondary instability, fortu- 
nately not too close to the threshold because a phase shift 
mode has not been disregarded [27 . In consequence, a 
model like this could be tried out for a Floquet expo- 
nent q = 2/3. Meanwhile, in the secondary bifurcation 
to ALT we show that Xalt/^s ^ 2, for this reason this 
instability belongs to the spatial period-doubling oscilla- 
tory bifurcation predicted by the CI theory, although it 



is build in the framework of a weakly nonlinear regime. 



B. Stability diagram 

Experimental studies for silicone oils with Pr > 
0(10) |40] show that each regime in the stability dia- 
gram is well defined because it is confined in large re- 
gions of the parameter space. In other words, there is 
a long pathway towards turbulence with a great diver- 
sity of regimes. Therefore consecutive instabilities can 
be evaluated separately. Our system {Fi^ 75) exhibits 
about eight different regimes in the explored space of the 
control parameters ((i, ATy) including bistable regimes. 

For d < 4 mm the quantitative study of the evolution 
of the amplitudes (l^sl and |A^±|) does not allow to dis- 
cover the nature of the transition from the ST pattern 
to the mixed ST/ ALT pattern . We may assert that for 
d < 4 mm the basic pattern ST never looses its stability. 
On the other hand, we have measured that the surface 
occupied by the localized domains in the ALT pattern 
with respect to the total surface of the spatiotemporal 
diagram is variable as we increase ATy. 

The system is bistable in the range of intermediate 
layers (4 mm < d < 4.5 mm), where a ST/ ALT do- 
main may coexist with a TW domain. Bistability is the 
reason why, in the secondary bifurcation to TW, two 
states with high spatiotemporal coherence are solution 
of the system for the same range of the control param- 
eters {d^ATy). In Fig. [tJ) it has been sketched how the 
mixed ST/ ALT pattern can be considered as a "transi- 
tional" pattern selected by the system for the subcritical 
bifurcation to TW. As we approach the bistable regime 
and until the system has fully transitioned to TW, the 
emerging traveling modes with wavenumbers k2s/3 play 
a destabilizing role towards the nearby oscillating modes 
with wavenumbers As the frequencies in the sub- 

critical interval are roughly the same for both oscillatory 
modes (see Fig. [sJd), cOsr^ close to the subcritical param- 
eter —Ec the stationary mode, Ms(/cs,0), activates the 
resonant ALT modes, M±y(kg^2i^sr ), before the TW 
modes, M±y{k2s/3,^sr)^ which are responsible for the 
next global bifurcation. Attending to the first nonlinear 
resonant criterion, a resonant interaction cannot be trig- 
gered between the stationary and the traveling modes. 
For the propagative modes in this mixed ST/ ALT pat- 
tern, we have measured (Fig. ^) a finite jump of wave- 
lengths {5k ■ d ^ 0.25) in comparison with the average 
wavelength of the propagative modes in the ALT pattern, 
this difference is smaller than the experimental error. 

As we get closer to the codimension-2 point the sub- 
critical parameter Ec tends to zero. At this critical point 
and depending on the way we approach it, the system 
may bifurcate to TW from two different patterns: ho- 
mogeneous (PC) or cellular (ST). From d ~ 4.5 mm the 
threshold for the PC towards the ST pattern gradually 
decreases until (i ~ 3, where the ST pattern settles with 
no threshold. This experimental result is in agreement 
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FIG. 12: (a) For = 4 mm, AT^ = 16 K) the spatiotemporal diagram of a bistable pattern in the secondary bifurcation to 
TW where the previous pattern ST/ ALT and TW coexist. The boundary domain where the stationary mode Ms is coupled 
with the traveling modes Mv± has been sketched (with /i = 0.37). Lsub is the average width of the boundary domain. The 
evolution of the amplitudes of the fundamental modes for the same spatiotemporal diagram (b) at t = 131 s and (c) t = 489 s. 
The growth of As corresponds to the resonant interaction at the front. 



with the theoretical prediction (a linear stability analy- 
sis) for a similar experimental setup [4r that shows that, 
as far as the Gaussian profile of the temperature field flat- 
tens a smaller threshold for the ST pattern is expected. 

From the PC pattern and for d > 4.5 mm, the system 
undergoes a supercritical bifurcation towards TW. At the 
threshold, we may think of a whole array of hotspots 
that cooperatively decide to travel towards a privileged 
direction, so the correlation length is supposed to diverge 
and consequently, for this supercritical transition, there 
are no fronts (even though the presence of fronts is not 
necessarily linked to subcriticality). 

From the results obtained for the secondary bifurcation 
to ALT, we can assert that bist ability is the cause for the 
coexistence of the patterns TW and ALT. Furthermore, it 
is the cause for the existence of localised drifting domains 
in the ALT regime in transients. This is the case for an 
abrupt upward jump (in terms of ATy) above the TW 
threshold where the presence of drifting domains in ALT 
have locally broken the symmetry imposed by the basic 
TW pattern (Fig.[5|. 

It should be pointed out that the ALT pattern 
is strongly nonlinear and far from the threshold it 
shows a typical turbulent dynamics of defects (see 
Fig. 11 ). In the threshold to the ALT pattern and above. 



the system always selects the two counterpropagative 
modes My±{ks/2^^ALT) which restore the stationary one 
Ms{ks^O) by nonlinear resonance. By contrast, at the 
threshold of the secondary bifurcation to TW there are 
two possible modes that may become unstable. My- or 
M^+, only one of them will be selected by the system due 
to nonlinear effects. 

In both secondary bifurcations the fronts that connect 
bistable patterns are stationary. Certainly, for intermedi- 
ate layers, the predominant thermograviational effect to- 
gether with the underlying cellular pattern (for ST/ ALT 
and ALT) are the major causes for the strong locking 
that restrains uniformly the diffusion of a new pattern 
into the original one at the asymptotic states. On the 



contrary, for thinner layers {d < 4 mm) the thermocap- 
illary effects probably account for the fluctuating fronts 
in the ST/ ALT pattern. 



C. Convective/absolute instability towards 
traveling waves and the subsequent alternating 
pattern 

For an experimental ID convective system which un- 
dergoes a secondary bifurcation to TW, we have shown 
how the existence of hysteresis can be tested studying the 
amplitude of the fundamental traveling modes (|^v+| or 
\Ay_ I ). According to the existence of a subcritical inter- 
val for the secondary bifurcation to TW, at (i = 4 mm the 
finite jump of wavenumbers and frequencies is produced 
at AT^^ = 17.2 K (Fig. [§i,b), before the jump of the 
amplitude of the mode My_ at AT^^ = 17.7 K (Fig. [T^) 
for an ascending sequence (with step 0.5 K). Close to 
the threshold these results might not be surprising, if we 
follow the usual form [42], under the symmetries of the 
problem (parity-breaking bifurcation), for the coupling 
between the slowly varying phase (j) and the amplitude A 
of the traveling mode: 



At — Ax 



. + /3A- 



(1) 
(2) 



where a and P are the coupling parameters. In our sub- 
critical bifurcation to TW, f{A) is expressed in terms of 
a fifth-degree polynomial of A. In accordance to equation 
Eq. ([T]), when we mechanically perturb the fluid the exci- 
tation of the amplitude of the traveling mode (Fig.|9| im- 
plies a coupling between the critical mode associated to 
the secondary bifurcation (M^+ or My-) with the phase 

gradient of the basic state ST, that is: At = a(l)xA-\ . 

The attenuation of the traveling mode in Fig. ^ shows 
that the domain in TW disappears because of a relax- 
ation dynamics for the unstable mode. 
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We have observed the bistable pattern (coexistence of 
the ST/ ALT and TW patterns) for positive and smah 
values of e because our measurements are restricted by 
the size of the step AT^. From the experimental results 
at the threshold, if we compare two ascending sequences 
with different steps 0.5 K (Fig. [t]) and 0.3 K, we find 
that the discontinuous jump in the amplitude \Ay_\ is 
produced before with the smaller step at AT^^ = 16.3 K 
{/ST^c = 17.7 K with step 0.5 K). This fact may suggest 
that for smaller steps we would measure decreasing values 
for AT^^, nevertheless we have shown that the subcritical 
interval increases the smaller the depth is. 

At d = 4 mm, for the two subsequent secondary bi- 
furcations the fronts, in the bistable regimes, are normal 
fronts that connect two stable states. We have shown 
that these fronts are stationary for a large range of the 
control parameter ATy. As ATy is increased the new 
unstable state will finally invade the whole cell, conse- 
quently the system has developed a global instability. 
This kind of behavior has also been found in catalytic 
reaction systems where stationary fronts appear in ID 
heating with a platinum wire [43]. Theoretically the 
minimum velocity of the front is achieved at the Maxwell 
point [29 , so in accordance to our results we are sup- 
posed to have a ''Maxwell intervar . This interval, with 
a minimum and constant velocity, has already been in- 
cluded in some theoretical models for reaction-diffusive 
systems [lH [45] . These models (for heterogeneous catal- 
ysis) show how a system with two competing diffusivities 
{D2/D1 <^ 1) is able to sustain a quasi-stationary front 
for a finite interval of the control parameter with veloc- 
ity Vp ~ \/^. According to the high Prandtl number 
in our experiment, if we apply the previous theoretical 
model, the thermal diffusivity would play the role of D2 
and therefore, for intermediate layers this model fits the 
stationary front dynamics. 

Because of a small linear group velocity, that has been 
measured below the threshold to TW for a FKPP front 
(Fig. [9|, the convective nature of the TW instability 
has an absolute character, therefore the threshold of the 
global instability is considered the same as for the con- 
vective one. Also, this global-convective/absolute insta- 
bility to the TW pattern compares well with the the- 
oretical work for oscillatory instabilities in the case of 
large Pr number fiuids in which the buoyancy effect is 
becoming important [46 . This is the case for d = 4 mm 
where Bod ^ 1-7. On the other hand, the subcritical 
behavior of the system vanishes for d > 4.5 mm, so from 
the codimension-2 point {Bod ^ 2.38) and for increas- 



ing depths no fronts do exist. For a supercritical model it 
has been proved that the "convective" critical control pa- 
rameter for a confined system verifies = + ^{^x'^)i 
where £a is the "absolute" critical control parameter for 
the corresponding extended system [47 , so for d > 4.5 
mm the approximation £c ^ is valid. 

The bifurcation to ALT has also a subcritical behavior 
and at the threshold the pattern ALT becomes unsta- 
ble at any position in the cell. The existence of a sta- 
tionary front could be attributed to an inhomogeneous 
heating along HL. Nevertheless, the subcritical interval 
is about 6 K and the inhomogeneity of the temperature 
field between opposite sides of Lx is less than 1 K. More- 
over the system bifurcates supercritically for d > 4.5 mm. 
The strong coupling between the different Fourier modes 
present for large values of ATy does not allow to charac- 
terize the successive bifurcations quantitatively. There- 
fore the quantitative analysis is useful in the weak non- 
linear regimes, at least for the secondary bifurcation to 
TW. 

For fiuids with high Prandtl numbers these results con- 
tribute to broaden the experimental evidence about cel- 
lular patterns which undergo subcritical transitions, so 
modelling systems with similar behavior could be carried 
out in the future. Close to the threshold of the secondary 
bifurcation to TW, the global unstable mode by nonlinear 
coupling drives the dynamics through a bistable regime. 
The domain that has become unstable towards TW has 
previously crossed a mixed pattern where the resonant 
triad has been activated. Far from the secondary insta- 
bility to the ALT pattern strong nonlinear interactions 
between fundamental and harmonic modes are the rea- 
son for a defect dynamics. This defect mediated chaos is 
the last stage of our bifurcation scenario, restricted by 
the available control parameter ATy. 
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